module geom
    integer :: nx, ny, nz, nxy
!    integer, pointer, dimension(:) :: nxs, nxe ! starting and ending x-coordinates of the core region for each y-direction
!    integer, pointer, dimension(:,:) :: nodel ! radial node index corresponding to the x- and y- cordinate
    integer, parameter :: GDW=1, GDE=2, GDN=3, GDS=4, GDT=5, GDB=6
    integer, parameter :: G11=1, G12=3, G21=2, G22=4, G1=1, G2=2
    integer, dimension(2,2) :: SCGG
    data SCGG /1,2,3,4/
    integer, dimension(4)   :: DRCOV
    data DRCOV /GDN,GDS,GDW,GDE/

contains
subroutine initgeom
    implicit none
    include "pspec.FI"
    integer :: i, j, l
    nx = ncmx
    ny = ncmy
    nz = 1
    nxy = nx*ny
!    allocate(nxs(ny))
!    allocate(nxe(ny))
!    allocate(nodel(0:nx+1,0:ny+1))
    
    ! assign proper geometric values
    l = 1
!    do j=1, ny
!        nxs(j)=1
!        nxe(j)=nx
!        do i=1, nx
!            nodel(i,j)=l
!            l=l+1
!        enddo
!    enddo
!    nodel(0,:)=0;nodel(:,0)=0;nodel(nx+1,:)=0;nodel(:,ny+1)=0
    
endsubroutine

endmodule
